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We investigate intersubband polaritons formed in the asymmetric quantum well (AQW) embed- 
ded into the semiconductor microcavity and study the eflFects of spin-orbit interaction (SOI) acting 
on intersubband excitations. The spin-orbit interaction of Rashba and Dresselhaus type remove the 
spin degeneracy of electrons with finite value of in-plane momentum and allow four types of inter- 
subband excitations. While optical spin-flip transitions are suppressed, the spectrum of elementary 
excitations shows the appearance of upper, lower and middle polaritonic branches based on spin- 
conserving transitions. The accounting of finite photon momentum leads to non-zero average spin 
projection of electronic ensemble in the first excited subband under cw excitation for both isotropic 
(Rashba) and anisotropic (Rashba and Dresselhaus) SOI. We predict the possibility of spin current 
generation in the considered systems with long coherence length. 



I. INTRODUCTION 

Intersubbband transition in semiconductor quantum 
well (QW) plays significant role in the modern optoelec- 
tronics due to numerous possible applications in optical 
devices operating in the infra-red and terahertz frequency 
domainSfi"— The dependence of energy distance between 
subbands on QW width allows to adjust the frequency 
of photon emitter or detector in relatively easy way com- 
paring to the usual interband transition. Moreover, the 
implementation of multiple QW samples gives the possi- 
bility to create devices with high efficiency, in particular 
quantum cascade lasers<^ 

An important characteristic of intersubband transition 
as compare to interband transitions is the peculiar opti- 
cal selection rules which only allow absorption of TM po- 
larized electromagnetic mode because dipole element of 
transition for TE polarized mode is zeroi^ Furthermore, 
it is possible to improve the efficiency of light interac- 
tion with absorbing media by placing it into the semi- 
conductor microcavity. This allows to achieve the strong 
coupling regime when in case of intersubband transition 
cavity photons are constantly absorbed and emitted and 
mixed light-matter modes are formed^"— 

Up to now, a consideration of the spin properties for 
intersubband polaritons was never performed. On the 
other hand, the spin electronics or spintronics is one of 
the most fastly developing areas of mesoscopic physics. 
The main issue of spintronics is spin transfer in the sys- 
tem and generation of spin currents. Being widely stud- 
ied nowadays, it proposes the devices which operate on 
other principles comparing to usual electronics jiiii^ for 
instance - the spin field effect transistor 

One should note, that contrary to the case of the inter- 
subband polaritons, the spin properties of interband cav- 
ity polaritons were widely studiedJ^ Moreover, by anal- 
ogy with spintronics, its optical counterpart, spinoptron- 
ics, became recently an area of the intensive studies.— 
In this domain the role of Rashba SOI is played by so- 
called TE-TM splittingi^ and optical analogs of various 



spintronic components were theoretically proposed and 
experimentally realizedJ^ 

One of the basic concepts in spintronics is the spin- 
orbit interaction (SOI )^^'^^ which appears in semicon- 
ductors due to intrinsic bulk inversion asymmetry (BIA) 
or structure inversion asymmetry (SIA). 

The part of SOI appearing in the systems with struc- 
ture inversion asymmetry, e.g. asymmetric quantum 
wells (AQW), is known as Rashba term. It can be repre- 
sented by introduction of the in-plane effective magnetic 
field causing the precession of the electron spini^ The 
corresponding Hamiltonian reads 



HsiA = a{axky - ayk^) = - {^siA ■ c), 



(1) 



where a is a Rashba SOI constant and ^sia = 
2ahr^{ky] —kx) denotes the effective magnetic field (mea- 
sured in frequency units). The diagonalization of the 
electron Hamiltonian accounting for Rashba SOI gives 
modifies dispersions where two different spin components 
arc split in energy for non-zero values of electron in-plane 
momenta. This removes the spin degeneracy in the sys- 
tem and makes possible effective spin control due to the 
possibility to tune the Rashba coupling parameter a by 
external gate voltage applied perpendicular to structures 
interfacei^"— 

Spin-orbit interaction arising from the bulk inversion 
asymmetry (BIA) is known as the Dresselhaus term. 
Similarly to the Rashba SOI, it leads to the appearance 
of linear in k effective magnetic field oriented in plane of 
the QW, but has different symmetry. The corresponding 
Hamiltonian can be written as^^ 



Hbia 



P{<7xkx - CTyky) = - {^BIA ■ Cr), 



(2) 



where /? is Dresselhaus constant for material, ^bia — 
2j3h~^{kx]—ky). In realistic QWs usually both types 
of spin splitting are present. This leads to strongly 
anisotropic pattern of the effective magnetic field acting 
on electron spin in the reciprocal space. 
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In this article we study the effect of Rashba and Dres- 
selhaus SOI on the intersubband optical transitions and 
formation of intersubband polaritonic states. Electrons 
in energy subbands of asymmetric QW are subjected to 
SIA and BIA spin-orbit interaction and both fundamen- 
tal and upper subbands are spin split for fc 7^ 0. This 
opens four different optical transitions and allows to form 
five different polaritonics states. In the article we show 
that optical transitions with spin-flip are suppressed in 
semiconductor microcavity and spin-conserving excita- 
tions interacting with cavity photon give birth to three 
strongly coupled polaritonic modes with peculiar spin po- 
larization. Tuning the pump conditions one can generate 
the spin currents with long coherence length. 

The article is organized as follows. In Sec. Ill Al we 
present a Hamiltonian for intersubband excitation sub- 
jected to Rashba and Dresselhaus spin-orbit interaction. 
In Sec. IIIBI the possibility of optical spin orientation 
due to non-zero photon momentum is discussed. In Sec. 
Ill CI we discuss spin currents generation by linearly po- 
larized light for intersubband transition. In Sec. IIIII we 
introduce the strong light-matter coupling in the system 
and show the spectrum of excitation with corresponding 
spin polarization of polaritonic modes. Finally, Sec. IIVI 
summarizes the results of the article. 



II. PHOTOABSORPTION OF INDIVIDUAL QW 
WITH SOI, SPIN POLARIZATION AND SPIN 
CURRENTS 

A. Model Hamiltonian 

1. Rashba SOI 
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FIG. 1: (Color online) (a) Sketch of the geometry of the sys- 
tem. The set of asymmetric QWs is bounded into cavity 
formed due to the effect of total internal reflection (waveguide 
geometry), (b) The schematic representation of four possi- 
ble transitions between upper and lower subbands splitted by 
SOI. A denotes bare transition energy, k+ and fc_ show the 
radii of Fermi surface for different SOI subbands. 



We consider a system of asymmetric GaAs/AlGaAs 
quantum wells embedded into the microcavity, where TM 
polarized cavity photons are confined (Fig. [D^a)). TE po- 
larized mode can be excluded from consideration, as it is 
not coupled to intersubband transition in the dipole ap- 
proximation. As compared to the case of symmetric rect- 
angular quantum well, the asymmetry introduces Rashba 
spin-orbit interaction which induces the spin flips for the 
electrons moving with finite value of the momentum k. 

The generic Hamiltonian for the considered system in 
secondary quantization representation can be written as 

(3) 

+ H.c] + ^ Efbl\ + ^ {gqai^^^,ak+q,2,ab\ + H.c), 

q k.q.cr 

where we have chosen the axis perpendicular to the in- 
terface of the QW as spin quantization axis z. a^jg-, 
ak,3,(y are creation and annihilation operators for electron 
with wave vector k and spin a in lower [j — 1) or upper 
(j = 2) subband, gq is the electron-photon interaction 



constant which originates from dipole matrix element of 
intersubband transition and can be calculated as^S 

9 = ! (4) 

cav) + Q 

where Lcav is cavity length, A is separation energy be- 
tween levels, uJo{q) denotes cavity mode dispersion, eo 
and e are vacuum permittivity and relative material di- 
electric constant, respectively, dio stands for the dipole 
matrix element of the transition and A is an area of the 
sample. 

The first term in Hamiltonian ([3]) describes free par- 
ticles and the second term corresponds to the Rashba 
spin-orbit interaction, where aj are Rashba coefficients 
for different subbands. j = 1,2. Note, that in general 
ai 7^ q;2.— This term contributes to the mixing of f and 
states. The third term is free cavity photons energy and 
fourth one describes the interaction between photons and 
intersubband transition which conserves electrons spin. 

It is convenient to diagonalize the electronic part of the 
Hamiltonian ^ introducing the new operators of spin 
states a = +, " oriented perpendicular to the direction 
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of the momentum k 



1 



Ck,± 



(5) 



^/2 

where Ok — arctan(fcj,/fca.) is the angle between radial 
vector k and x axis. The Hamiltonian of the light-matter 
coupling written in new operators reads 

k,j,(T q k,q 

X 4 1 +Cfc+q,2,+ + .9 — {k,q)cl 1 _Cfe_q,2 - + .g+-(fc,g)x 

(6) 

^ 4,i.+c/c+q,2 - + .g-+(A:, q)cl .^ _Ck+q,2.+]bl + H.C., 
where Ej {k) stands for the standard Rashba dispersions 



E,{k)=E,{k)±af\k\ 



(7) 



shown schematically in Fig. [TJb) for upper and lower 
subband j. 

The electron-photon coupling coefficients now arc func- 
tions of electron momentum k and are spin-dependent: 



g++ik,q) = gqc- 



9+-{k,q) = igqe 



sm 



'+{k,q) = igqe 



g — {k,q) = gqe 
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9qg^ 



gqg^ 



(8) 



(9) 

- .99.9-+, 
(10) 

- gqg — , 

(11) 

where gq is an absolute value of light-matter coupling 
coefficient defined above. 

Note, that for non-zero photon momentum 9 7^ the 

spin- flip transitions become possible , g \. 7^ 0) 

and all four transitions with different energies are al- 
lowed (see Fig. (Ub)). However, the main impact for 
photoabsorption comes from electrons situated close to 
the Fermi surface, for which 9 ^ kp. Hence spin-flip 
matrix elements are proportional to the momenta ratio 

g_l , g |_ ~ q/kp ~ 0, their oscillator strength is small 

and only two transitions H > + and >■ — can be ex- 
pected to be visible in the experiment. Note, that the 
energies of these transitions are different, and thus they 
can be excited selectively by tuning the frequency of the 
excitation beam. 



asymmetry the Drcsselhaus SOI term is usually present 
in the Hamiltonian. The combination of both Rashba 
and Dresselhaus SOI leads to spin anisotropy and pecu- 
liar spin orientation in the subbands. The total electronic 
Hamiltonian H = Hsia + Hbia (see Eqs. H]) and ([2])) 
can be diagonalized in the same fashion as above by in- 
troduction of the operators 



Cfe,± 



1 

71 



±- 



^a^ + /32 + 2a/3 sin 0k 



.(12) 



The dispersion relations for electrons arc not cylindrically 
symmetric and have non-trivial form 



Ej{k) = Ej{k) ± \k\^a] + 13] + 2a/3sin26'fc, (13) 

where j denotes the number of subband. 

In full analogy with the case when only Rashba SOI 
is present, one can get the matrix elements for all four 
optical transitions. While the calculation is straightfor- 
ward, their explicit expressions are rather cumbersome 
and we do not present them here. Similarly to the case 
of isotropic Rashba SOI the spin-flip transitions are sup- 
pressed, w 0. 



B. Spin polarization 

1. Rashba SOI 

The achievement of non-zero spin polarization is one of 
the main goals of spintronics. Regarding the optical gen- 
eration of spin polarization (spin orientation), there are 
many proposals based on the interband excitation, in par- 
ticular by circularly polarized light 3^ As well, there are 
several proposals for spin generation with intersubband 
transition, where optical selection rules imply that only 
linearly polarized light is absorbed and achieving of non- 
zero spin polarization becomes a formidable task. One 
mechanism of generation is based on different strength 

of H > + and — — > — transitions due to valence band 

mixing which modifies transition matrix elementsi^ An- 
other approach requires different SOI effective field in 
both fundamental and excited subband leading to over- 
all non-zero spin polarizatioui^ Here we want to study 
the possibility of spin polarization generation which is 
linearly proportional to the photon wave vector 9. 

The average spin of the electron in the second quanti- 
zation representation reads 



-,29.30 



'k.(T 



Ck., 



(14) 



2. Rashba and Dresselhaus SOI 

In the realistic QWs Rashba SOI is not the only in- 
teraction acting on the electron spin since due to bulk 



where we defined the spin operator S = ^{crx,cry,crz) 
with (Ti being Pauli matrices given for each direction (i = 
x,y^ z). It is reminiscent to the first quantization defini- 
tion of average spin given as Sk.± = (4'fc.± I'^'l^'fe.i)'^ 
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First, let us consider the case of fully thermalized elec- 
tron gas and calculate the average spin within subband. 
The average spin projection onto x and y directions for 
the case of Rashba SOI are given by 



(+1^.1 

{+\Sy\ 

{-\s.\ 

(-\SJ 



^ Jo (27r)^ Jo 

= ^(-COS^^,)= / / d0,(-COS0fe), 

^ Jo (271-)^ Jo 

= ^(-sin0,)=/ tS/ defc(-sin0,), 
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^ kdk f^"" 

= } cos9k= / t;— / d9kCOs9k, 
V -'o (271-)^ Jo 



and one can see the orientation for spin in H — h and 

transitions given by vectors S+ — (sin^fe, — cosfl^, 0) and 
S_ = (— sin^fc, cos 0) with three spatial components 
S = {Sx, Sy, Sz)- The modified radii of Fermi surface for 
+ and — spin subbands are written as k± = fcF(l T ?7-r) 
(sec Fig. mb)), where rji^ = -pfj- The electron spin 
orientation for spin subbands is shown in Fig. ^a). Ob- 
viously, the integration over angle 9k yields zero average 
spin projections, {zL\Sx,y\i) = 0. 

Now, let us consider the situation when electrons are 
constantly optically excited to the upper subband due 
to cavity photon absorption. The spin of an electron in 
the presence of spin-orbital interaction of Rashba type 
is perpendicular to its momentum direction. There- 
fore, with accounting of finite wave vector of the photon 
one can see that spin state in Fig. [2{c) is no longer 
given by vector S , but S. This corresponds to the 
change of spin states to S+ = (sin 9i;+q , — cos 9k+q , 0) and 
S_ = (— sin 9 k+q, cos 9 k+q,0). Now the boundaries of in- 
tegration are changed due to shift of Fermi circles and for 

selective excitation of H — I- and transitions are given 

by expression 



Y sin^ i9k 



k^ cos{9k - 9q) 



(15) 

The average spin projections with accounting of pho- 
ton momentum can be rewritten using relations 

k cos 9k + qcos9q 

- 2kqcos{9k — 9q) ' 



cos ( 



sin 9 k+q 



k sin 9k + q sin 9q 
y/k^ + q^ + 2kqcos{9k~eq)' 



(16) 



where 9q denotes angle between vector q and x axis, and 
all variables are shown explicitly on the sketch in Fig. 
[5Jc) . One should keep in mind that Fermi circle for q 7^ 
is 0fc-dependent. Therefore, it requires advanced momen- 
tum integration with consequent angular integration and 
can be done only numerically. To obtain the dimension- 
less quantity connected to average spin of electron gas 
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FIG. 2: (Color online) (a) Orientation of electron spin in 
-I- and — Rashba splitted subbands without accounting of 
photon momentum (g = 0). (b) Vector plot of spin orientation 
with accounting of finite photon momentum (q/kp = 0.3, 
9q = 0). (c) Schematic representation of non-zero photon 
momentum q influence on spin polarization. 



(or spin polarization) one should normalize the result of 
integration dividing it by population in each subband 



n± 



1± 

ilT ■ 



One sees that while for g = case average spin is 
zero, the accounting of Fermi circle shift gives prefer- 
able direction to the spin of excited electrons which is 
[010] for 9q = incident angle shown in Fig. [IJb). This 
effect is similar to spin-galvanic effect observed previ- 
ously for interband case and connected to non-zero charge 
current >i2i^ 



2. Rashba and Dresselhaus SOI 



Following the scheme used in the previous subsection 
for the spin-isotropic Rashba splitting, we derive the 
same quantities for spin-anisotropic case of combined 
Rashba and Dresselhaus interaction. 

The spin orientation can be found using the spin states 
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for Rashba and Dresselhaus case (|12p and is given by 

a sin Ok + (3 cos 9^ 
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where we defined dimensionless Rashba and Dresselhaus 
constants a = and = ^j^-g, respectively. The 

vector plot of spin orientation for different Rashba and 
Dresselhaus constants is shown in Fig. [3] (a, b). One sees 
that while for a = 2/3, (3 = 1/3 the spin pattern is close 
to Rashba case but with preferred directions being [110] 
and [TTO], for equal Rashba and Dresselhaus SOI strength 
{a ~ P ^ 1/2) the spin orientation approaches step-like 
function which changes sign in 0^ = 37r/4 and 9k = 77r/4 
points. The Fermi surfaces for upper (+) and lower (— ) 
spin subbands now are given by expression k± = fci?(l 4= 

r^{9k)) with r^{9k) = + + 2al3 sm29k. 

Accounting for non-zero photon momentum result into 
change of spin states S(6'k) — > S (^k+q)- Additionally, 
one should account for angle dependent shifted Fermi 
surfaces 




(c) 



FIG. 3: (Color online) Orientation of electron spin in -I- and 
— Rashba and Dresselhaus splitted subbands for equal SOI 
strength (a) and for 5 = 2/3, /3 = 1/3 case (b). (c) The 
spin vector plot for equal Rashba and Dresselhaus SOI with 
accounting of finite photon momentum which leads to Fermi 
surface shift in [110] direction. The areas / (yellow filling) 
and // (pink filling) show the difference in spin orientation. 



k'^ = kF{l+v{9k^q)) Jl- 



with r]{Qk,q) 
Ok,q = arctan 



cos 



-dq) 



(17) 



sm Uk —Q sm t 



:^a2 -1-/32 + 2q;/3 sin 28^,9 and 



Similarly to the isotropic 



-q cos t 



k cos 6 

Rashba case, for q = the average spin vanishes, while 
for finite photon wave vector q the net spin is not equal 
to zero. The effect is most clear for equal Rashba and 
Dresselhaus SOI. Then the average spin can be seen as 
difference of areas / and // in Fig. [31[c). Moreover, one 
can note that 5 = /3 = 1/2 case is strongly anisotropic re- 
garding the photon incidence direction 9q (Fig. b)). 
While for 9q = n/A it has maximum absolute value, spin 
polarization does not appear in the system for 9q = Sn/A 
angle of incidence. The average spin vector (S) is ori- 
ented perpendicular to q. Dashed lines in Fig. |4l^a, b) 
correspond to not equal Rashba and Dresselhaus SOI and 
depict transient regime between spin-anisotropic to spin- 
isotropic case. 

We calculate the dependence of spin polarization on 
the orientation of an incident light angle 9q. The cor- 
responding polar plots are represented in Fig. for se- 
lective excitation of (a, green solid circle) and -|--|- 



(b, purple solid circle) transitions. Note, that in case 
where only Rashba term is present, the absolute value 
of average spin is the same for all 9q. The orientation 
of average spin vector is defined by the 9q angle and is 
aligned perpendicular to the vector q (Fig. 21 (S) label 
on the plot). 

Furthermore, we study the dependence of spin polar- 
ization on the wave vector ratio q/kp which defines the 
relative value of the Fermi circle shift under optical exci- 
tation. Since photon wave vector value usually does not 
overcome inverse micrometer q < 10~® m~^, the exper- 
imental way to tune q/kp ratio is changing the carrier 
density. The plot of net spin polarization as a function 

of momenta ratio q/kp is given in Fig. Sl^c) for and 

++ transitions (green and purple solid lines). Higher 
values of spin polarization for lower spin subband come 
from the fact that A:_ < kp < k^. One can see the linear 
dependence for the case of isotropic SOI for the realistic 
values of wave vector ratio. In the inset of Fig. |4l^c) we 
show the behavior of average spin on the large scale. It 
shows saturation for region where photon wave vector ap- 
proaches the value of Fermi circle radius with spin polar- 
ization approaching unity. However, firstly, this situation 
does not have substantial physical background since for 
small concentrations {q kp 10~^ m~^) even small 
excitation leads to the strongly non-equilibrium situa- 
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FIG. 4; (Color online) (a) The polar plot of absolute value 
of an average spin vector as a function of photon incident 
angle 9q for -I— I- transition. Green solid circle corresponds 
to the Rashba SOI only case (5 = 1, /3 = 0). Blue solid 
circles show strong anisotropic dependence of spin polariza- 
tion on incident angle 6q for equal Rashba and Dresselhaus 
strength (5 = /3 = 1/2). The blue dashed line shows tran- 
sient behaviour for 5 = 2/3 and /? = 1/3 case, (b) The polar 
plot of absolute value of average spin vector as a function of 

photon incident angle 9q for the lower spin transition. 

Purple solid circle corresponds to the isotopic Rashba SOI 
case (5 = 1, /3 = 0). Red solid circles show spin polariza- 
tion as a function of incident angle 6q for equal Rashba and 
Dresselhaus strength {a = /3 — 1/2). The red dashed line 
shows transient behaviour for 5 = 2/3 and P — 1/S case. All 
dependences are calculated for q/kp ~ 0.06. (c) The spin po- 
larization dependence on the wave vector ratio q/kp for the 
Rashba case (green and purple solid lines) and Rashba with 
Dresselhaus case (blue and red dashed lines) for 5 = 2/3 and 
P = 1/3. The inset shows the behaviour of spin polarization 
for Rashba only case on large scale of argument. 



tion and this situation can not be treated within our 
approach. Secondly, the visible separation of polaritonic 
modes implies that separation between bare transition 
modes Cs = 2(a2 — ai)kp is greater than photon decay 
rate (F sa 1 meV), which puts the lower boundary for 
kp value. Finally, for the ratio q/kp closer to unity one 
should account for H — and — h transitions on the same 
footing with -|--|- and transitions, and we do not con- 
sider this situation. 



C. Spin currents 

In the following section wc briefly discuss the possi- 
bility of spin current generation by intersubband optical 
excitation. The concept of spin currents lies in the origin 
of spintronics since it allows to design devices based on 
manipulation of spin degree of freedom of the carriers. 
However, effective creation of spin currents is not a triv- 
ial task. The proposed schemes include injection from 
ferromagnetic leads to semiconductor creation of non- 
equilibrium distribution of carriers by applying voltagoi^ 
or by optical excitation^Sri^ in structures with spin-orbit 
interaction. The photocurrents associated with non-zero 
spin polarization were studied in the context of interband 
transition with spin-split bands shined by circularly po- 
larized lightj^i^i^ The inverse effect namely the appear- 
ance of charge current due to non-equilibrium population 
of electron spins was observed by Ganichev et al^ and 
named spin-galvanic effect. The review on the subject 
can be found in Ref. [s^. 

In the system with spin-split upper and lower subbands 
the intersubband transitions occur with spin conserva- 
tion. Due to the difference in the SOI constants for upper 
and lower subband, the H — h and transitions have dif- 
ferent energies. An absorption of TM polarized photon 
leads to the electron in the upper subband with momen- 
tum k (first we consider vertical transition) and "hole" 
in the lower subband with the opposite momentum (see 
Fig. [S]). Moreover, we account for the Rashba SOI with 
peculiar spin orientation in the spin subbands. The ex- 
citation of carries on the Fermi disc leads to the genera- 
tion of spin and charge currents (Fig. [51 top). However, 
for the case of vertical transition charge current is zero 
and situation corresponds to the generation of pure spin 
currents. The expression for the pure spin current in the 
general form can be described by the pscudotcnsor J with 
components 



„36 



k 



TpTr 



^vX(k)' 



(18) 



where v(k) is velocity operator and p(k) is spin density 
matrix, is relaxation rate. This formula applied to the 
intersubband transition describes the spin currents in v 
direction with spin aligned in /i direction. 

In the previous section we have shown that accounting 
for the finite photon momentum leads to the non-zero 
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FIG. 5: (Color online) Sketch of spin current orientation (j, 

black arrows) for the (a) and ++ (b) transition. The 

coloured lines show the spin orientation of electrons (s). 



average spin of intcrsubband excitation. Thus, it creates 
a non-equilibrium situation similar to the electric field 
displacement of Fermi disks when the charge current 
in the system does not vanish and spin current with spin 
polarization defined by incident angle 6q is generated. 
While we do not derive the expressions for spin current in 
q 7^ case, the behavior can be qualitatively understood 
in connection with previous section. 

Finally, the implementation of strong-coupling be- 
tween photonic mode and intersubband excitation affects 
the spin currents. Due to coherence of cavity mode the 
resulted spin currents will have the long coherence length 
and can be used for spinoptronic applications 4^ 



III. INTERSUBBAND POLARITONS WITH SOI 
A. Elementary excitation spectrum 

In the case when QW is placed inside a microcavity 
which is tuned in resonance with intersubband transi- 
tion, the photons can undergo multiple re-emissions and 
re-absorptions and hybrid intersubband polariton modes 
can be formedi^s^i The powerful theoretical tool for their 
description is many body diagrammatic techniqucJ^ The 
quantities which define the physical properties of the sys- 
tem are photon Green function and polarization opera- 
tor, shown diagrammatically in Fig. ^a,) and Fig. EJb), 
respectively. The poles of the renormalized photon Green 
function give the spectrum of elementary excitations in 
the system, while imaginary part of polarization operator 
allows to calculate the absorption spectrum. 

If one neglects the electron-electron interactions, the 
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FIG. 6: Diagrammatic representation of intersubband polari- 
tons. (a). The Dyson equation for Green function of microcav- 
ity photon interacting with H — h and intersubband tran- 
sitions II±,±. g± corresponds to the interaction constant, in- 
dices 1 and denote the excited and fundamental subband, 
respectively, (b) Dyson equation for polarization operators 
governing the absorbtion spectrum of intersubband quasipar- 
ticles coupled to the cavity mode where He = 11-1-+ + 11 



polarization operators corresponding for 4 possible types 
of transitions in the system ,-| — , — h) read 



TTO-CT 

^^0 



(2) 

k-t-q.cr' 



(19) 



Z7 



where we defined ip = 9k — 0k+q and a, a — ± denote dif- 
ferent spin subbands (-1- and — ). The energy of electron 

in fundamental subband is E^\k) = ft^fc^/2m ± Q;i|fc| 
and in upper subband is (fc+g) = A+h'^{k+q)'^ /2m± 
a2\k + q\ and dimensionless coefficients g~~' (g, ip) are de- 
fined by Eqs. ©-HI]) and read 



\9++iq,^)\ ^ l.9l_(9,¥')l -cos^ v/2= 7^(1 + cos 

(20) 



\9-+iq,v)\ 



2 
1 



sin = -(1 — cos(^). 



It was already mentioned in Sec. Ill Al that while \g±±\ 
interaction constants corresponding to H — h, tran- 
sitions are close to unity, the -I — , — I- transitions do 
not play substantial role. The quantitative estimate 
can be done using realistic values of electron concentra- 
tion with Fermi momentum of the order kp sa 10^ m~^ 
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(n = 10 cm~ ) and microcavity photon momentum 
q in the point of anticrossing q = 10^ m~^ for nega- 
tive detuning (5 = — 5 meV. The angle is typically small 
If = Ok — 0k+q ~ a.rcsm{q/kF) ~ 10~^ which reduces the 
osciUator strength of spin-flip transitions gi^, (g, ip) by 
five orders of magnitude. Consequently, they give minor 
contribution in the case of strong light-matter coupling, 
and we omit them in the Dyson equation for cavity pho- 
ton Green function shown in Fig. |6|a). The polariza- 
tion operators corresponding to H — h and transitions 

can be calculated by performing analytical integration on 
angle in (jl9p with subsequent numerical integration on 
absolute value of wave vector k. 

The dispersion of elementary particles in the system is 
defined by the poles of renormalized photon Green func- 
tion D{q,Lu), which can be found by solving the Dyson 
equation, and reads 



D 



Dn 



1 



(21) 



where Tl±.± arc polarization operators described before, 
qq is given by Eq. Q and Dq is a bare photon Green 
function, 



-Do(w,(7) 



2hu>Q{q) 



h'^uj^-h^uji{q) + 2iToJo{q)' 



(22) 



where uJoil) denotes cavity mode dispersion and F is a 
broadening of photonic mode due to the finite lifetime of 
the cavity photon (taken to be « 10 ps). The explicit 
expression for the electron-photon interaction matrix el- 
ement gq is given by Eq. (|4]). 

The spectrum of the elementary excitations of inter- 
subband transitions coupled to the cavity mode with ac- 
counting of Rashba SOI is plotted in Fig. [TJa). Here 
we considered AQW of L = 12.8 nm width with fre- 
quency of bare transition equal to A 100 meV at /c = 
point and concentration of electron gas n = 10^^ cm^^. 
The Rashba constants for lower and upper subband were 
taken to be ai ~ 0.9 mcV nm and a2 = 6 meV nm, 
respectively)^^ and we neglected Dresselhaus terms for 

simplicity. Both -| — h and spin-dependent transitions 

are in strong-coupling regime with the cavity mode. As 
transition frequencies are different, one observes the for- 
mation of three polariton branches. This differs from the 
usually considered spin-independent case for which only 
two polaritonic branches exist. 

The energy splitting between ++ and transitions 

is approximately eg = 2{a2 — ai)kp. For instance, with 
the given parameters a2 and ai and for concentration 
n = 10^^ cm^^ it has value eg = 0.8 meV, which is 
already comparable to the broadening of the photonic 
mode. Therefore, for the experimental observation of 
the effect one should search for the sample with large 
difference in Rashba constants, and InAs based QWs are 
promising candidates for this purpose. 

While the consideration of simultaneous Rashba and 
Dresselhaus interaction can be done in the same fashion, 
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FIG. 7: (Color online) Dispersions of intersubband polariton 
modes with accounting of isotropic Rashba splitting. The 
red line corresponds to lower polariton (LP) and coincides 

with excitation dispersion for large momenta, while blue 

line shows upper polariton (UP) which is reminiscent to -I— |- 
excitation for small wave vectors. The green line describes 
middle or mixed polariton (MP) being a transitional branch 
between lower and upper polaritonic branch. 
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FIG. 8: (Color online) The spin polarization degree as a 
function of pumping laser energy for the photon wave vec- 
tor q — 6.25 jj.m^^ . Changing the laser frequency allows one 
to excite selectively upper, middle or lower polariton modes. 
This results into change of the net spin polarization. 



it does not lead to qualitative differences comparing to 
the isotropic Rashba case and we do not address this 
point here. 

Intersubband polariton modes can be selectively ex- 
cited by the resonant excitation of a given frequency. As 
it was discussed above, such excitation will induce spin 
currents in the system. The net spin polarization as a 
function of the excitation frequency is shown in Fig. [S] 

Due to the coupling with photonic mode, spin currents 
induced in the coupled QW-microcavity system will have 
much bigger decoherence length compare to the electron 
currents 4^ This makes the considered system a promising 
candidate for spintronic and spinoptronic applications. 
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IV. CONCLUSIONS 

In conclusion, we analyzed the optical properties of 
spin-dependent intersubband transition in asymmetric 
quantum wells. We have shown that accounting of fi- 
nite photon momentum leads to the optical orientation 
effect. We also calculated the spectrum of elementary 
excitations arising from strong coupling of the photonic 
mode with an intersubband transition of an asymmetric 
QW. The calculated spectrum of elementary excitations 



show the appearance of upper, lower, and middle polari- 
tonic modes. The possibility of the generation of spin 
currents with long coherence lengths is discussed. 
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